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ABSTRACT 

Preston-tube data have been obtained on a sharp-nose, ten-degree cone In 
the NASA Ames 11-ft TWT and in flight tests. During analyses of the laminar- 
boundary-layer data , errors were discovered In both the wind-tunnel and the 
flight data. The apparent errors in the 11-ft TWT data are relatively minor 
and were easily corrected. However, the errors in the flight data are much 
more severe. A great deal of effort was expended in the search for a rational 
procedure for correcting this data. A correction procedure is recommended 
which forces the flight data to exhibit some of the orderly characteristics of 
the wind-tunnel data. 

Subsequent to correcting the wind-tunnel data, a correlation is developed 
between Preston-tube pressures and the corresponding values of theoretical 
laminar skin friction. Because of the uncertainty in correcting the flight 
data, a correlation for the unmodified data is developed, and, in addition, 
three other correlations are developed based on different correction procedures. 
Each of these correlations are used in conjunction with the wind-tunnel 
correlation to define "effective" freestream unit Reynolds numbers for the 
11-ft TWT over a Mach nimber range of 0.30 to 0.95. Using the preferred 
correlation, based on the recommended rearrangement of the flight data, the 
maximum effective Reynolds numbers are approximately 6.5% higher than the 
normal values over a unit Reynolds number range of 9.8 to 16.4 million per 
meter. These maximum values occur between freestream Mach numbers of 0.60 
and 0.80. Smaller values are found outside this Mach number range. These 
results indicate wind-tunnel noise J^fects the average laminar skin friction 
much less than it ^fects boundary- layer transition. 

Data on the onset, extent, and end of boundary-layer- transition are 
also summarized for these tests. The wind-tunnel data indicate a Reynolds 
number, based on distance to end-of- transition, is a unique function of noise 
and Mach number. 

Finally, a procedure is described for studying the relative effects of 
varying nose radius on a ten-degree cone at supercritical speeds. Preliminary 
results indicate increasing nose radius promotes boundary- layer transition 
and separation of laminar boundary layers. 
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C^. Skin-friction coefficient (2 

C Nondimens ional difference between theoretical and correlated 

' skin^-friction coefficient - ^f,c^^^f,t^ 

C Pressure coefficient based on the difference between a Pitot 

P and static pressure reading [{Pp - P^)/q„,] 

(Cp)^ Pressure coefficient on surface of cone [(P^^ - P»)/0 
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Nondimensional effective height of Preston tube (2Yg^^/h) 

L Axial length of cone, 44.5 in. 
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Mach number 

[ 

Friction Mach number, {\/yP^^)^ 

Static pressure at outer edge of boundary layer 

Preston tube pressure 

Static pressure at wall 

Difference in pressure between a Preston tube and wall pressure 
Freestream dynamic pressure (pU„^/2) 

Reynolds number based on Uq and external diameter of a circular 
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Reynolds number based on diameter of a sphere and freestream 
conditions 

Reynolds number based on surface length from stagnation point 
and V evaluated at the reference temperature of Sommer and 
Short (UgS/V) 
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Critical Reynolds number for a sphere 

Reynolds number based on external diameter of a circular Preston 
tube and properties at outer edge of boundary layer, U^d/v^ 
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X 
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X * 
AX 


Reynolds number for incompressible flow thru a pipe of 
diameter D, (Uj^O/v) ' 

Freestream unit Reynolds number per meter (U^v^) 

Effective freestream unit Reynolds number per meter 

layer^and^r*^^*" ^^sed on properties at outer edge of boundary 

Reynolds number based on freestream properties and X^ 

Surface distance measured from a stagnation point 
Temperature 


Velocity parallel-to bounding surface 

Nondimensional velocity used in the law-of-the-wall (u/U ) 

Velocity parallel to axis of cone 

Velocity at outer edge of boundary layer 

Mean or average velocity in a pipe flow 

Velocity calculated from Preston-tube data and P 
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Classical wall -shear-stress velocity (lyp)**'* 

Freestream velocity 

oa^fn^l ® flattened Pitot probe in a directio 

parallel to the wall but normal tot he undisturbed streamlines 

Distance along axis of cone 

difference for incompressible flow 

logioLAPp d^ZA pv^] 

mlasuJpd of cone; Appendix B, distance along axis 

measured from apex of a sharp-nose cone 

Most forward station at which Preston-tube measurements began 
fayer"transition”'"'^^‘^® boundary- 

UyertranJition^''"^®^^^ boundary- 


Dimensionless pressure difference for compressible 
batic flow log,o(Up 

Allen's correlation parameter logio(L) d/v') 

P 

Length of boundary-1 ayer-transit ion zone 


nonadia- 


xi 


y Distance measured normal to the wall; Appendix B, denotes 

perpendicular distance from axis of cone 

y,. Distance of yeonietric center of Preston-tube from wall 

y' Nondimensional distance from the wall as used in the law-of- 

t he -wall (U^.y/o) 

y* Dimensionless shear stress for incompressible, isothermal 

flow logio[r d*/4pv^ or 0.25(U d/v)^J 

W 1 

Y Effective height of face of Preston tube = height of an undis- 

turbed streamline above the wall which-has_a_ total pressure 
equal to the measured Pitot pressure 

Y* Dimensionless shear stress for compressible, nonadiabatic 

flow iogi.(y*ff/p„V> 

Y.* Allen's correlation parameter loguC*^ U d/v') 


Greek Letters 

a Angle-of-attack, defined to be positive for nose up 
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6 Boundary layer thickness 

/ 

6* Displacement thickness for compressible boundary layer 

0 Momentum thickness for compressible boundary layer 

)i Molecular viscosity 

V Kinematic viscosity 

p Density of fluid 

t Shear stress at wall 
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w ' w w ^ 
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I. INTRODUCTION 


A. Background 

Since the 1920's, it has been known that freestream turbulence within 
wind tunnels caused differences between ostensibly identical tests in dif- 
ferent wind tunnels and in flight. It was discovered, during these early 
years, that these differences in low-speed wind-tunnel tests could be ac- 
counted for by defining an "effective" freestream unit Reynolds numbe^L 
which was higher than the corresponding flight unit Reynolds number. The 
now classical procedure for defining an "effective" freestream Reynold num- 
ber was developed based on the drag of a sphere. This procedure is based 
on the well-known fact (even in the 1920's) that the dra^ coefficient of a 
sphere (and a circular cylinder) drops sharply at a critical Reynolds num- 
ber of the order of 10®. This sharp drop in drag coefficient is caused by 
transition from a laminar to a turbulent boundary layer near the shoulder 
of the sphere. The turbulent boundary layer remains attached overa longer 
distance al-ong the surface of the sphere and thereby reduces form drag. 

According to Pope and Harper^ a convention was eventually adopted to 
define the critical Reynolds number as the one which corresponds to a sphere 
drag coefficient of 0.30, which occurs roughly at the maximum in d(Cg)/d(Ro) . 
Careful measurements in free-air showed that the critical Reynolds number was 
approximately 385,000. The values measured in wind tunnels were always found 
to be less. This led to the formulation of a turbulence factor defined as 
the ratio of the free-air critical Reynolds number divided by the correspond- 
ing wind-tunnel value, i.e., 

TF = 385,000/(Re^)y.^ . 

The "effective" tunnel unit Reynolds number is then defined by 

^*^®ft^eff ' ^^®m ^ • 

Additional details of this procedure can be found in Ref. 1. Unfortunately, 
this method becomes inaccurate above = 0.35 as compressibility effects 
become important. In particular, the sphere generates a rather strong shock 
slightly downstream of the shoulder of the sphere at high subsonic Mach num- 
bers. This phenomena changes the relation between drag and Reynolds number, 
e.g. , see Lu’’ . 
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In low-speod tunnnls 0.50) the freestreani unsteadiness 1s pre- 

dominantly vorticity which Is controlled in modern tunnels by a carefully 
desiqned turbulence-management section upstream of the test section, e.g., 
Loohrko and Nagib* and Eckert, et al.** Whereas, in transonic tunnels the 
doiiiin.int source of flow unsteadiness is wall -generated noise caused by flow 
through and across the ventilated test-section walls which are required to 
establish steady_f3ow near Mach cme. The-noisc intensity levels in transonic 
tunnels typical peak between = 0.70 and 0.80, e.g., Reed, et al.'’ The 
AEDC Boundary Layer Transition Cone (see Fig. 1) was designed to calibrate 
the effects of noise on boundary-layer transition. It is superior to a 
sphere because a sharp-nose cone does not generate a significant transonic 
shock.* The need for such a calibration device was indicated by discrepancies 
between numerous transonic wind-tunnel tests of models at ostensibly identi- 
cal flow conditions. A particularly wel 1 -documented study of differences in 
static aerodynamic data has been obtained with the same model of a C-5A trans- 
port aircraft in three major transonic wind tunnels; the results have been 
reported by Treon, et al.^ The differences between the three different sets 
of tunnel data were reduced by accounting for "relative" Reynolds number ef- 
fects between facilities. The AEDC 10-deg cone was used to define the dif- 
ferences in "relative" Reynolds number. As noted by Dougherty and Steinle: 

"These results substantiated the need for developing a method for 
predicting these corrections to Reynolds number to improve the 
extrapolation of wind-tunnel test results to full-scale flight 
conditions, i.e., a "turbulence factor" for transonic tunnels. 

Although the "effect" associated with these differences in transi- 
tion characteristics between tunnels is of prime importance in ad- 
justing the data, the "cause" is of particular significance since 
it relates directly to predicting the "effect" in new and different 
facilities. This illustration of improvement in agreement of re- 
sults between transonic facilities suggests the use of transition 
Reynolds number for such corrections to be both technically appro- 
priate and productive." 

In 1971, an extensive test program was begun in which the AEDC-BLT cone 
was tested with the same probe-traversing mechanism and instrumentation, see 
Fig. 1. (Additional details of the cone are given in Refs. 7 and 8.) In 
addition, the cone was accompanied by Dougherty from AEDC in order to assure 
tost procedures were matched as closely as possible. The primary purpose of 
this program was to simply detect the location of the boundary-1 ayer-transi- 
tioM-?one on the rone in different wind tunnels and in free-flight but at 

*Se'o T^jTpetiBT x' IT ToV^fiTrth'er d i scussTon dT^'TFu^Tf fects of nose bluntness on 
tratisonic flow about cones. 



Figure 1. AEDC Boundary Layer Transition Cone 


idonticnl (nr compnrablo) values of ROm and The technique of using 
a traversing Pitot probe to detect boundary layer transition has been used 
since the 1930's. In fact, the first Wright Brothers' Lecture by Jones^ in 
1937 describes the utility of this technique in flight tests. Figure 2 is 
taken from this paper and shows clearly the change in total pressure across 
the transition zone along different paths parallel to the bounding surface. 

The results of tests in six different tunnels is reported in Ref. 7, and 
a summary of results obtained in twenty three tunnels and a concluding flight 
test are reported by Dougherty and Fisher.® In this concluding report, 
Dougherty and Fisher found that the data for transition Reynolds number, 
based on the product of local unit Reynolds number and distance from nose to 
end-of-transition (X^), was proportional to with an_error band of 

+ 207.. This represents a significant step forwaxd_in the development of a 
procedure to calibrate flow qualUy in transonic, wind tunnels. 

The purpose of the work reported herein ts- to investigate the possibi- 
lity of using the traversing-Pitot-probe data to derive more precise and/or 
additional information concerning the effects of noise on flow quality. The 
basic approach, which was selected to achieve this objective, is to interpret 
the surface Pitot-probe data as Preston-tube data, i.e.., total pressures near 
the wall which can be related to skin friction. Unfortunately, four different - 
probes were used during the series of wind-tunnel tests, and two were used 
during the flight tests. Replacement was necessary because of probe wear 
(along the underside of the probe where it made contact with the cone), 
damage and/or deterioration during use. This did not introduce any signifi- 
cant problem with regard to detection of transition but does become very im- 
portant if the data are to be interpreted as Preston-tube data. In order to 
set the stage for this type of analysis of the data, the basics of Preston 
tub(?s and their use to measure skin friction is now introduced. 


TRANSITION BEGINS 
DISTANCE ALONG SURFACE 


VARIATIONS OF 

TOTAL PRESSURE THROUGH A REGION OF TRANSITION 
FROM LAMINAR TO TURBULENT BOUNDARY-LAYER FLOW 



^ stpjl Jil-bpA 

According to Proston^®, the British engineers Stephens and Haslani^^ 
suggested in 1938 that it should be possible to use the data from a Pitot 
tube traversed along a surface to infer skin friction. Apparently, this 
idea was not pursued until Preston's work during the early 1950' s. He 
developed a correlation between skin friction and the total pressure as 
measured with circular pitot tubes resting on the inside wall of. a pipe. 

In order to develop this correlation, Preston assumed the classical 
law-of-the wall is valid across the face of the probe and chose the- 
characteristic length to be the height of the geometric center of the 
probe above the wall, i.e., d/2. As shown in Ref. 12, this leads to the 
following relation between Preston-tube pressure and skin friction. 


<'‘p - 
pV2 


G[ 


^w(d/2V 

PV2 


or inversely, 


~^w d^ _ pr APp d^ i 

4pv2 ~ '' '■4pv2 •' 


( 1 ) 


This relation suggest a convenient method to determine skin friction .since 
the shear stress is uniquely related (for given fluid properties p and v) 
to the difference in pressures between a Pitot tube (OD = d) and a static- 
pressure orifice at the wall. Using Eq. (1) as a guide, Preston obtained 
measurements inside a pipe flow with circular Pitot tubes having four dif- 
ferent external diameters but constant ratios of Internal to external 
diameters of 0.6. Pipe Reynolds number was varied over the range lO"* < 

RCp <10^, Skin friction was determined via measurements of pressure drop 

over a known length of constant diameter pipe, viz., t = (Pi - P2)D/4L. 

w 

An empirical fit of the data led to the following correlation. 

y*^-1.396 + g- X* [Preston, 1 954] (2) 

Where y* 1og,o C^wdVAnvM and x* = logic (APpdV4pv^ , 

In this same paper, Preston also reported measurements with the tubes on the 
floor of d wind tunnel in slightly favorable and strong adverse pressure 
gradients. The use of Eq. (2) to convert the pressure measurements into 
skin friction coefficients led to a set of data that was consistent with 




6 


corresponding boundary- layer surveys with a ininature, flattened Pitot 
probe. 

Thus, in addition to the surface Pitot measurements with circular 
probes, Preston also used a small {h = 0.015 cm.) flattened Pitot probe 
to measure velocity profiles. In an effort to obtain accurate profiles, 
he employed an earlier estimate by Young and Maas^^ that the effective 
center of the probe would be displaced 0.25h toward the region of higher 
\telocity. This correction is ostensibly to account fora velocity gradient 
across the height of the probe face. This is needed because as noted by Chue^** 

"Errors caused by the presence of a shear flow across the mouth 
of the pitot tube are due to the following two effects: 

(i) the stagnation pressure is proportional to the 
square of the velocity and when this is integrated 
over the orifice, it will have a higher value than 
the stagnation pressure calculated from the square 
of the velocity at the geometric centre of the ori- 
fice; and 

(ii) the presence of the probe in a velocity gradient 
causes deflection of the stream lines toward the 
region of lower velocity. This deflection causes 
the probe to indicate an impact pressure in excess 
of that existing at the same location in the absence 
of the probe. 

The existence of the second effect has been qualitatively demon- 
strated by smoke photographs. Both of these effects are therefore 
seen to cause the probe to read high, which explains the outward 
displacement of the effective from the geometric centre." 

The net effect is the measured pressure corresponds to the total pressure 
of a streamline which is above the geometric center of the probe. Although 
Preston attempted to correct for this effect of shear, he noted that there 
is an additional effect associated with proximity to the wall. In particular, 
when the flattened probe was within a distance 3h of the wall, he found the 
displacement of the effective center was reduced and appeared to be a func- 
tion of Uxh/v. 

Preston did not attempt to define the functional relationship between 
y^, h and Reynolds number. However, this was subsequently undertaken by 
MacMillan'® for circular Pitot probes. His measurements indicated that 
when the geometric center of a circular probe is more than two diameters 
away from the wall, the effects of the wall on displacement of the effective 
center is nil, and the displacement due to shear alone is 0.15d (i.e., 

= 0.65d), independent of Reynolds number. He also established that the 
effect of the wall is to reduce this displacement and move the effective 
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probe position closer to the wall. This is easily seen when one realizes 
that the downward flow across the face of the probe (caused by shear) 
is iiiij)eded as a probe approaches the wall. Furtherinore, the oncoming 
streamlines (below the center of the probe) begin to lift upward and move 
over and around the probe instead of passing underneath between the probe 
and the wall, MacMillan proposed a single curve for a velocity correction 
which is to be added to -the measured velocity in order to account for wall 
displacement effects. The correctionJs a function only of y^/d and is 1.55K 
of the measured velocity when y^/d »0,5 and is zero when y^/d >2.0. This 
correction for wall effects is to be added to the displacement effects of 
shear. By expressing the wall effect in terms of a fraction of the measured 
velocity, MacMillan was able to define a correction which is independent of 
Reynolds number. However, since the measured velocity at a given value of 
y^^/d is a function of Re^, the displacement caused by wall proximity is also 
a function of Re^. Thus, MacMillan concluded that the total displacement of 
the effective center is a function of y^/d and U^d/v when 0.5 < y /d < 2.0. 

It is relevant to here note that any displacement effects, which occurred in 
Preston's data, are buried within the empirical coefficients of Eq. (2). 

In 1964, Patel published the results of an extensive set of tests 
with fourteen different circular Pitot probes and three different pipe 
diameters. He obtained a more accurate calibration for Preston tubes and 
established limits on the pressure-gradient conditions within which his cali- 
bration can be used with prescribed accuracy. Patel obtained empirical equa- 
tions for y* = f(x*) over three ranges of y*: (1 ) 3. 5 <y* < 5.3, (2) 1.5<y*< 
3.5, and (3) y*<1.5. These three regions correspond, respectively, to the 
fully-turbulent, the buffer or transition zone, and the viscous-sublayer 
regions of the classical law-of-the-wall . In incompressible flow, the normal 
Preston-tube Reynolds number range corresponds to the buffer zone, and for 
this region Patel obtained 

y* = 0.8287 - 0.1381 x* + 0.1437(x*)2 - 0.0060(x*)^ (3) 

where: 1.5 < y* < 3.5 or 5.6 < U.jd/2v < 55. Patel claims this correlates 

his data to within f 1.5% of x^. 

In the viscous-sublayer region, Patel found his data was correlated by 

y* = 0.5x* + 0.037, ( 4 ) 

when: y* < 1.5 or U,d/2v <5.6 . 
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In this near-wall region, the classical law-of-the-wul 1 exhibits the linear 
relation 

:: u/Ut ^ U-jy/v s . (5) 

In order to relate Eqs. (4) and (5), Patel introduced and defined the 
"effective" center of a round Pitot tube to be at 


. ( 6 ) 

By definition of the effective center, the velocity recorded by a Preston 
tube, Up, is the true velocity at 


1 


1 


- - iPp - j- - 

If this is substituted into Eq. (5), the results are 


p 2 


y* = 0.5x* - 0.5 logic (0.5 


(7) 

( 8 ) 


Now equating Eqs. (4) and (8) and solving for a value of 1.3 is 

obtained. 

Patel noted that this agrees "precisely" with the value determined by 
MacMillan. However, this value for is equivalent to MacMillan's re- 
sults for displacement due to shear alone, i.e., when y^/d > 2.0. Whereas, 
In the case of Patel's use of Pitot tubes resting on the bounding surface, 
the value of y^/d is 0.5. Thus, Patel's results for appears to be for- 
tuitous. but an argument can be made that makes this plausible. Firstly, 
MacMillan noted an additional displacement correction for viscous effects 
on Pitot probes in zero-shear leads to larger, positive displacements when 
Ufd/v < 25.** Secondly, as previously discussed, wall proximity effects 
result in a reduced or a negative displacement. Therefore, it appears 
(assuming no significant experimental errors) that Patel's value for 
occurs because viscosity and wall effects cancel each other. 


**As IfToVed in 'the re“v^w article by'^^e"^. It is generally agreed that the 
pressure coefficient (Cp) for a Pitot probe is greater than one when U d/2v 
300. However, there is no concensus between existing experimental data and 
theoretical results as to precisely how Cp varies with Reynolds number. 
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II. RESULTS OF PRELIMINARY ANALYSES 
OF WIND-TUNNEL DATA 

The traver«i1ng Pitot probes, used during the AEOC-BLT Cone tests, 
are of the flattened or ovaUshaped type. Figure 3 shows typical dimen- 
sions of the probe used during flight tests of the cone. The probes 
used in wind-tunnel tests are geometrically similar^. The particular 
probe, used during tests in the NASA Ames 11 -ft Transonic Wind Tunnel 
(TWT), had a height of 0.0097 in. 

Since Patel's results are for circular Preston tubes, they cannot 
be applied directly to the AEDC Cone tests. In addition, these tests 
were conducted at transonic speeds and compressibility effects are ex^ 
pected. With regard to flattened Preston tubes, Quarmby and Das^’ con- 
ducted an experimental study and calibration of six oval.-shaped probes 
when used as Preston tubes. When x* > 4.6, they found these probes gave 
exactly the same calibration relation between y* and x* as was obta-ined 
by Patel (Eq. 3) if the external height of the probe face is used in place 
of d. At lower values of x*, the negative displacement of the effective 
center caused by wall proximity was larger (=5%) for the flattened probes 
with aspect ratios between 1.5 and 1.9.** The following calibration equa- 
tion correlated the measurements of Quarmby and Das to within 1.5% of t . 

w 

y* = 0.5152 + 0.1693 x* + 0.0651 (x*)^ 

for 3.38 < X* < 6 (9) 

Since these results for oval -shaped Preston tubes agree so closely with 
Patel's results and Patel's value for =1.3 appeared to be appropriate 

in the visco us-sublayer of a turbulent wall-flow, it was initially decided to 
use this same value in an attempt to correlate the traversing Pitot probe 
data obtained within the laminar boundary layer on the AEDC Cone. This ap- 
peared to be reasonable in light of the fact that the x*'s for the cone 
data were > 5.5. Although this is equivalent to assuming is independent 
of Mach number, Reynolds number, velocity gradient across the face, and 
aspect ratio, this assumption was attractive because it greatly simplified 
the analytical work. 

Now turning our attention to compressibility and Mach number effects, 

**^This is consistent with the idea that fTbw about the face becomes more 
two-dimensional as aspect ratio increases and more of the flow passes up 
and over the face rather than around the sides. 
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DIMENSIONS OF TRAVERSING PRESTON-TUBE PROBE USED DURING FLIGHT TESTS 




AMon'" has perfonnod the most comprehensive analysis of Pres^-on tubes In 
supersonic boundary layers. He developed a correlation using three inde-- 
pendent sets of simultaneous measurements of Preston "tube pressures and 
skin friction via a floating-element force balance. These data were ob- 
tained within flat-plate, turbulent boundary layers and with Mach numbers 
in the range: 1.6 < < 4.6. Allen selected the same basic dimension- 

less parameters as Patel; except he chose to evaluate the fluid properties 

p and V at a refecence temperature developed by Sommer and Short^®, and 

the velocity Up was calculated from Pp and tha wall pressure (=» Pg) using 
standard compressible flow relations.** 

xj 5 log,, I ^ } - log,, (Up d/v') (10a) 

5 log,, I Rej(p'Cp/Pp)^|. log,, (.7 U^d/v') (10b) 


The primes denote properties evaluated at the Sommer and Short reference 
temperature, viz. , 


1 ? 


T7T^ = 0.55 + 0.035 
e e 

The correlation derived by Allen is 

+ 0.45 yr^ . 

(11) 

Y* = -0.4723 + 0.7814 X* 

+-0.01239(X*)“ . 

(12) 


[Allen, 1977] 


Allen found that the majority of the skin-friction-coefficient data were within 
+1 5rt to -12% of Eq. (12). This rather large scatter, compared to the in- 
compressible pipe-flow calibrations of Patel and Quarmby and Das^’, is 
at least partly associated with the much greater sensitivity and vulnerability 
of floating-element balances to extranneous errors. tt 

Obviously, the parameters used by Allen are logical candidates in 
any attempt tc correlate the transonic cone data. However, the basic pur- 
pose of a reference temperature is to permit use of skin friction formulas 
for incompressible flow to estimate compressible skin friction by evaluating 
fluid properties at the reference temperature. Thus, the resulting reference 
properties represent an "average" value across a boundary layer. Whereas, 
small Preston tubes encounter only the flow near the wall. Therefore, it 
appeared to us that properties based simply on the wall temperature would 
be more apropo. The utility of evaluating properties at both of these tempera- 


**The details can be found in the report by A11en^“ or Reed, et al.'* 
HAllen ‘ has discussed the various error sources in floating-element force 

ins^riment^^Ref^ zV suggested an improved design for this type of 


turos was invest1(jated» and the results are reported following a summary of 
the wind tunnel data. 

Wi nd Tunn el Data 

Although the AEDC-BLT Cone has-been tested in twenty three different 
wind tunnels, only the anaXyses of subsonic data from the NASA Ames 11 -ft 
Transonic Wind Xunnel (TWT) is reported herein. Suface Pitot-probe surveys 
were taken along the cone between axial stations 10 and 89 cm. aft of the 
nose. Table I lists the twenty one subsonic flow conditions at which the 
cone was tested.** The pattern of typical pressure surveys at high and low 
Reynolds numbers are shown^ respectively, in Figures 4 and 5. 

As a point of departure, this research seeks a correlation between the 
Preston-tube pressures, measured within the laminar portion of the cone's 
boundary layer, and the corresponding theoretical values of skin friction. 

If successful, the intent was, and is, to compare such a correlation with 
the corresponding correlation for flight data. It is thought that such a 
comparison can lead to the definition of an "effective" unit Reynolds 
number for the 11 -ft TWT. However, a literature search for Preston-tube 
data within laminar boundary layers turned up only one reference, viz., 
Prozorov.^^ 

Prozorov obtained surface Pitot-probe measurements within low-speed, 
flat-plate, laminar boundary layers. He used four circular Preston tubes 
and three rectangular-shaped probes with aspect ratios of 5.21, 5.21, and 
5.00. Although his data exhibited considerable scatter, he concluded, 
based on his measurements, that is a universal function of Upd/v (or Up 
h/v) for both laminar and turbulent boundary layers and is independent of 
probe geometry. Within the accuracy of his data, has nearly a constant 
value of 1.3 when R^ >100 and approaches 1.8 when R^j =32. This essential- 
ly verifies the results of Preston for turbulent pipe flows. However, 

Prozorov found T^d^/pv^ to be a different function of R^j for laminar boundary 
layers compared to what Preston found. The two calibration curves for 
T dVp^J^ diverge for Rj > 100. He explained this in terms of the following 
equation which was derived from a Maclaurin series expansion of Up at small 
distances from the wall and the conservation of mass and momentum for steady, 
two-dimensional flow. 

^Although surveys were also conducted at low supersonic speeds, the shock 
Interactions between the traversing-probe assembly and the hemispherical 
Pitot-static probe were not completely calibrated, and thus the supersonic 
data were not included in this research. 

■•■^his is inconsistent with the results of MacMillan'^ which indicated Keff 
for circular Preston tubes is a function of U-rd/v and the results of Quarmby 
and Das'’ which indicated Keff is also a function of aspect ratio. 13 
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Prozorov claims this is valid in laminar, transitional and turbulent flows 
provided the probe height is within the viscous -sublayer. According to 
Eq. (13). \ = fn(Rjj) if and only if is also and the pressure gradient 

is negligible. _ _ 

It was decided not to use Prozorov 's equation (R^, dP^dx)] 

for the development of a correlation of laminar Preston-tube data for the 

following reasons, 

1. It was not expected that the_probes used in the subject tests 
would always be within the viscous sublayer. 

2. For efficiency, it is desirable to have a common procedure for 
analysis of both the laminar and turbulent data. Use of the 
classical law-of-the-wall leads to a procedure which is applicable 
to both types of flow, provided pressure gradients are small. 

3. The static pressure gradients along the cone surface were known 
to be small, and Patel's work^® defined the errors to be expected 
when using-calibrations of Preston tubes, based on the law-of-the- 
wall, in turbulent flows with small pressure gradients. 

Thus, for analyses of the laminar data, it was decided to use the basic cor- 
relations parameters utilized by Patel and Allen for turbulent wall-flows. 


Computation of Boundary Layer and Da ta Analysis 

The distribution of static pressure along the surface of the sharp cone 
at subsonic speeds is assumed to be defined by the inviscid theory of Wu and 
Lock.^" Predictions for pressure coefficient along the surface of a 10-degree 
cone are shown in Figure 6 as a function of freestream Mach number. This in- 
formation and the known tunnel freestream conditions are used to calculate 
flow conditions along the outer edge of the boundary layer. The conical 
laminar boundary layer is then calculated using a computer program developed 
at Stanford University by Crawford and Kays^^ which they have labelled STAN-5. 
The resulting distributions of laminar skin friction and boundary layer pro- 
perties are then matched with the corresponding values of surface - Pitot 
measurements . 

It was arbitrarily decided to only use Preston-tube data at 1/2 in. 
intervals beginning with the most forward station at which data were obtained. 
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This resulted in a total of 148 data points along the cone at the various 
and listed in Table 1. The following linear equation was first used 
to correlate the Preston-tube measurements with the corresponding values of 


theoretical skin friction, 

Y* * A X* + B T* + C (14) 

where 

Y* = log.o 

X* 5 log.o (Up 05b) 

T* : log,o (T7T^) (15c) 

W " '^eff 


The reference temperature was introduced to account for small departures 
of the fluid properties p and v from the wall values. The coefficients A, 

B, and Cjwere determined by a least-squares fit of the data. This resulted 
in the following semi -empirical correlation. 

Y* = 0.655X* + 2.095T* - 0.895 (16) 

= Z* - 0.895 

A plot of Eq. (16) is presented in Fig. 7 along with the indiitidual data 
points. The corresponding differences in skin friction coefficient are shown 
in Figure 8. The rms value of = (C^ t “ ^f c^^^f t 6.7%. 

Next, an equation quadratic in X* was tried. The resulting correlation is 

Y* = 0.273(X*)2 - 2.618 X* + 1.645T* + 8.92 _ (17) 

A plot of this equation and the associated data are shown in Fig. 9. The 
corresponding C*^ is presented in Fig. 10. When the correlation parameters 
of Allen (Eq. 10) are used to fit the same data, the rms value of is 8.6%. 

Thus, the parameters defined in Eqs. (15) appear to be superior for correlating 
this particular data. An examination of Fig, 10 reveals that the data in 
the upper left-hand corner is separated from the majority of the data. These 
six data points corresponds to Run No. 57 (M^, = 0.80, Re„, =48x10®). When they are 
deleted, the rms value of is reduced from 6.7% to 5.2%. Finally, the seven 
data points, which form almost a vertical line on the left of Fig. 9, correspond 
to Run No. 72 (M^^ = 0.80, ROj^ =13x10®). When these points are also deleted, 
the correlation becomes 
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Y* = 0.0942(X*)=^ - 0.438X* + 2.023T* + 2.272 . (18) 

The associated ms value_of Is 4.93%. 

These results are good compared to the correlation of Allen^® but are 
rather. large compared to the very small scatter {= 1%) of the correlations for 
incompressible flow of Patel (Eq. 3) and Quarmby and Das (Eq. 9). Although 
greater scatter may be expected for compressible flows, somewhat less scatter 
is expected- for a correlatian_af the subsonic cone data because errors associ- 
ated with floating-element balances are not present as they are in the data 
considered by Allen. Thus, the question arises: how can the data be better 

correlated? This led to a reexamination of. the data and the development of 
an improved correlation- when is treated as a variab-le. The details of 
this second-analysis of the wind tunnel data are discussed in the next section. 


III. REANALYSIS OF WIND TUNNEL DATA 
USING VARIABLE 

Reexamination of the papers by McMillan,^® Patel, and Quarmby and 
Dasi7,26 igj yg conclude that, In general, the effective center of a 
Pitot tube is a function of: 

•<eff “ l<eff (^T-h/v, y^/h, w/h) . (19) 

In the case of a Preston tube, y^/h = 0.5, and aspect ratio (w/h) Is a 
constant for a given probe. When these restrictions apply, Eq. (19) re- 
duces to Kg^^(U^h/v). Since, In general, wall shear stress is a function 
of Reynolds number, pressure gradient, Mach number and heat transfer, we 
can expect for a given Preston tube to also be a function of these 
variables. If this conclusion Is true, it is necessary to Interpolate 
from the STAN-5 boundary- layer profiles. This has been done by finding 
the position within the theoretical laminar profiles at which the total 
pressure is equal to the measured Pitoto pressures. Table II provides a 
summary of the results for each wind-tunnel flow condition. 

It may be noticed that only nineteen cases appear in Table II as com- 
pared with twenty one in Table I, — In the process of tabulating (C )pj^g> it 
was discovered that the exact flow conditions for Run Nos. 27.411 and 39.545 
were in doubt. After a brief attempt to ascertain the correct values it 
was decided to drop these two cases -from further consideration. 

In addition to and Pp, Table II also includes noise measurements 
^^pVms obtained with a 0.635 cm microphone mounted flush with the 
surface of the cone at a distance 45.7 cm aft ofthe nose and 135 degrees 
around from the Preston tube. As discussed in the introduction, Dougherty 
and Fisher® have correlated boundary layer-transition with this type of noise 
data. Thus, it is anticipated that this data will be relevant to the defini- 
tion of an "effective" unit Reynolds number for the 11-ft TWT. 
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TABLE II (cONT'd) 

WIND TUNNF.L DAIA USED IN DEVELOPMENT OF I.AMINAR CORRELATION 





4.5 2115.3 0.977 

5.0 ?'’ 93.9 0.999 

5.5 2072.5 1.012 

5.0 < 053.0 1 . 0?3 



TABLE II (cONT'd) 

DEVELOPMENT OF LAMINAR CORRELATION 
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Xiie inothod used to define Kgfr^ has the effect of adjusting the height 
above the wall at which Pp is measured. This procedure is expected to lead 
to an improved correlation between Pp and because the Preston-tube pressure 
is forced to be consistent with the theoretical boundary-layer profile and skin 
friction. However, a high or low value of Pp and for a given value of 
and h leads to a numerically different relationship between X* and Y*. Higher 
values of Pp produce a more nonlinear correlation. 

In order to make this point apparent, it is helpful to see a representa- 
tive graph of the axisynimetri c laminar boundary layers in standard law-of-the- 
wall coordinates, i.e., u'*' vs. y"^. Three typical velocity profiles are shown 
in Fig. 11 for = 0.60 and three diffet::ent unit Reynolds numbers. The key 
observation to note from this figure is that Y^^^ “ 35,. and at this 

height the normalized velocity n* is only 7% below the linear relation u’*' = 
y'*’. When a Preston tube is completely submerged in this linear region, Patel'® 
has shown that a linear relation between y* and x* results, e.g., see Eq. (4). 
As the probe height increases and/or boundary layer thickness decreases, the 
relation between y* and x* becomes more nonlinear, e.g., see Eq. (3). In. the 
case of the subject data, a relatively small nonlinearity is expected. 

The above discussion of how a correlation is influenced by high or low 
values of Pp, naturally leads to the question of accuracy of the measured 
pressures, and how can erroneous data for a given wind-tunnel condition be 
identified? This can be qualitatively assessed by comparing the corresponding 
values of for a particular case with the distribution ot for the 
majority of the data,** For this purpose, has been plotted as a function 
of U^h/v^, and and is shown in Fig. 12. It is relevant to here note 
that the slightly favorable pressure gradients are negligible over the range 
0.09 < X/L < 0.26 for which laminar Preston-tube data are available, see Fig. 

6 and Table II. Thus, the systematic variations in are apparently caused 
by changes in flow about the face of the probe wiith changes in: (1) Reynolds 

number, (2) Mach number, and (3) tunnel freestream disturbance levels. -These 
variations in effective probe height must be properly accounted for if a single 
correlation equation, with constant coefficients, is to be uniformly ivalid 
with respect to Mach number. 

The variations of Keff, shown in Fig. 12, show that~the effective height 

of a probe decreases as lUh/v increases and/or Mach number increases. In 

^ w 

either case, the pressure difference APp ^ Pp - increases. Since U^h/Vyy 
increases as a given probe moves forward toward the nose, It is obvious that 
**Here we assume the bulk of tRF^ata provides a valid reference. 
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Figure 11 












AP will increase also.' In the case of Mach number it is not quite so 
obvious that APp increases with In order to illustrate this phenomena, 

Kic). 13 presents the effect of M^ on the distributions of total pressure 
across the calculated laminar boundary layers on a 10° cone. Here, we see 
that for a given wall shear stress the total pressure increases with . 

This is true for approximately the inner half of laminar boundary layers. 

Based on the theoretical boundary layer calculations and the given probe 
height (0.0246 cm) used during these wind tunnel tests, the range of the 
ratio of probe height to boundary layer thickness is 0.34 < h/6 < 0.77. 

However, the effective height varies in such a way that Yg^^/5 is always less 
than 0.5. In fact, for the three cases shown in Fig. 13, Yg^^/6 is < 0.32. 

The individual values of and APp are listed in Table III. 

The above described effects of Mach number have also been observed by 
Bradshaw and Unsworth^’ in their analyses of data for supersonic, flat-plate, 
turbulent boundary layers. In addition, the fact that the effective height 
of a probe is less than the actual height explains why Allen' s^® Preston- tube 
and turbulent skin friction data for supersonic turbulent boundary layers 
appeared to follow the logarithmic law-of-the-wall even when d/6 was as large 
as 0.70. 

Only two subsonic wind-tunnel conditions were repeated, viz., M^ = 0.7 
and 0.8 at a Reynolds number of 13 x 10®. Comparisons of for each of 
these cases indicate a difference of 0,075 for M„ = 0.7 and 0.15 for M^ = 0.80. 
These differences translate, respectively, to differences in measured pressure 
of 1.1 kPa (0.16 psi) and 2.3 kPa (0.33 psi). Since the full-scale range of 
the pressure transducer used in the probe is 34.5 kPad (5 psid), the corre- 
sponding percent erri-.^s in pressure are 3.2% and 6.6%, respectively. These 
values are a measure of the repeatability and precision of the Preston-tube 

data. 

Since the distribution of for a given M,^ is expected to be contin- 
uous, the discontinuities between the data for unit Reynolds numbers of 9.8 
and 13 are also a measure of precision. The 11-Ft TWT was shut down between 
the runs for different unit Reynolds number, and individual Mach number 
cases were run in the order listed in Table 1. However, there were two 
exceptions to this order. The tunnel was started for run numbers 44-47 and 
was shut down after wards. The second exception occurred for run numbers 

*0nTy^TaTa“fTo1ir>Vh n is being used in this work. 

The Preston-tube measurements of Prozorov^ ’ for incompressible flat-plate 
boundary layers also exhibit a decreasing effective probe height with 
increasing U^h/v^, when IJ^h/v^ is less than 100. 
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THEORETICAL DISTRIBUTION OF TOTAL PRESSURE ACROSS COMPRESSIBLE 
LAMINAR BOUNDARY LAYERS ON A SHARP NOSE TEN DEGREE CONE 



Figure 13 
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70 (M - 0.7, lUi - 13 X 10") nnd 72 (M,,. - 0.8, Ko - 13 x 10'’} which 
Wf-To imrlnniK'd rtl. a hi(|hor unit Reynold?! number iiiinied lately after the 
preceedino runs (f>fi-61) for a lower Reynolds number. Thus, for tliese 
two rims, It is suspected that the pressure transducer was beinq 
influenced by unsteady temperatures and may not have achieved an 
equilibrium temperature. This phenomenon may have also contributed 
to errors in pressure measurement for other cases. For example, the 
K for run number 44 (M = 0.95 Re == 16.4 x 10") appear to be low. 

Using the entire 136 values of shown in Fig. 12, a 

correlation between Preston-tube pressures and theoretical, laminar 
skin friction can be calculated. The equation obtained from a least- 
squares fit of a quadratic to the data is 

Y* = 0.0174(X*)=’ + 0.3274 X* - 0.0392 T* + 0.4333, 
for 5.4 < X* < 6.3 and 1.0. 


The rms error in C^; ^ is now 0.901. This amount of scatter is 
comparable to the pipe flow calibrations of Patel and Quarmby and 
nas’ 


However, the lack of continuity between the data for Re^^ ^ 9.8 x 10" 
and the rest of the data is of concern. A reexamination of the data 
sheets for this Reynolds number indicate some confusion as to the correct 
gain factor for the X-Y plotter. For example, on the Preston-tube data 
sheet for Run No. 56 - 0.90, Re^,^ = 9.8 x 10"), the listed gain factor 

is 0.289 psi/cm. Use of this gain factor results in a Preston-tube pressure 
higher than freestream total pressure! A review of the tunnel testing 
procedure led to the conclusion that Runs 56-61 should be shifted so that 
the distribution of for Run No. 58 (M _ = 0.70, Re^^ = 9.8 x 10 ) 

fonns a continuous curve with Run No. 21 (M^ = 0.70, Re^^ = 13 x 10 ). 

In addition. Runs 70 and 72 were shifted, as a group, so as to align Run No. 70 
with Run No. 21. A plot of the revised K^^^'s arc shown in Fig. 14‘. The 
corresponding correlation Is 


Y* = -0.0136 (X*)’ + 0.6977(X*) + 0.1051 T* + 0.6669, 
for 5.4 • X* • 6.3 and M 1.0. 


Till' curves, wliich appear in I ig. 14, will be discussed in the following 
section on analyses of tlight data. 
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Figure 14 
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This equation is shown in Fiq. 15 with the individual data points superimposed. 
The nils error in C. is now 0.97?'., and the individual values of are 

presented in Fig. 16. Equation (21) is considered to be our best estimate 
of a correct correlation for the 11-Ft TWT. 

It may be noted that the contribution of the T* term is very small because 
the temperature ratio t'/T^ is near one and the log of 1 is zero. This 
simply indicates negligible departures of density (n) and viscosity (y) from 
values based on wall temperature. In anticipation of the need to compare the 
wind-tunnel correlation with a correlation based on flight data, it is obvious 
that such a comparison will be easier if T* is dropped. Upon doing this and 
refitting the data, the following simplied correlation is obtained. 

Y* = -0.0103 (X*)2 + 0.6653 X* - 0.5946, ^22) 

for 5.7 < X* < 6.3 and < 1. 

The rms error in C- has increased only slightly to 0.98%. Equation (22) 

* * ^ ' L. 

is plotted on Fig. 17 and is compared with an analogous correlation for the 

unmodified wind tunnel data. It should be noted that the small increase 

in rms error of C. , relative to the unmodified data, is caused by the fact 
■ » c 

that a given difference in leads to a greater difference in Preston-tube 
pressure as the height above the wall is increased, recall Fig. 13. 

It is here emphasized that the numerical values of and the 

coefficients in Eqs. (21) and (22) are valid only for the Ames 11-Ft 
TWT and the particular probe used during these tests. The numbers are 
expected to be different for different wind-tunnel environments and for 
probes with significantly different aspect ratio and/or face geometry. 

In particular, the coefficients are believed to contain information on the 

freestream disturbance levels which are peculiar to the 11-Ft TWT. With 

the attainment of this objective, a-xorresponding analysis of the 
flight data was begun. This is reported in the next section. 
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IV. ANALYf.IS or FLIGHT DATA 


The AFOC-ULT Cone was mounted on the nose of a McDonnell-Doucjlas 
1-15 aircraft for flight testsV^*^ Although the cone apex was over 2m('“ 

7 Ft) ahead of the aircraft nose, it was recognized that the flow field 
about the aircraft could influence flow about the cone. Hence, it was 
necessary to obtain measurements of static pressure distribution along 
the cone by employing a facsimile cone. This cone has the same exterior 
dimensions as the BLT Cone but was constructed so that interchangeable 
inserts in the wall permitted measurements of either static pressures 
or wall temperatures along the surface of the cone. The inserts for 
static pressure measurements had an orifice with an inside diameter 
of 0.107 cm {0.042 in.) Table IV lists the locations of these orifices 

along the surface of the cone. 

Most of the flight data were obtained for nonzero pitch and/or yaw. 

Since we are unable to calculate non-ax i symmetric boundary layers with 
STAN-5, only the cases having negligible a and 3 were selected for 
development of a Preston-tube/skin friction correlation. A_total of 
nine cases were Identified as having both « and 6 of the order of 0.1 
degree or less. A summary of the selected data is tabulated in Table V. 

The numbering system used to designate the various flights is-as follows: 

XXX • 

Flight Time of 

No . Day 

The listed values of freestream Mach number, unit Reynolds number, dynamic 
pressure, and pitch and yaw angles are the time-averaged values based on the 
measured conditions which existed during the time the Preston tube was 
within the laminar boundary layers. 


TABLE IV 

PRESSURE INSTRUMENTATION ON FACSIMILE CONE 


(a) Surface Distribution of Static Orifices 


fice 

X(in.) 

xy.k. 

d)(DEG) 

I 

7 

.157 

180 

2 

16 

.360 

0 

3 

17 

.382 

0 

4 

19 

.427 

0 

5 

20 

.449 

0 

6 

21 

.472 

0 

7 

22 

.494 

0 

8 

25 

.562 

0 

9 

26 

.584 

0 

10 

27 

.607 

0 

n 

28 

.629 

0 

12 

29 

.652 

0 

13 

30 

.674 

0 

14 

30 

.674 

90 

15 

30 

.674 

180 

16 

30 

.674 

270 

17 

31 

.697 

0 

18 

32 

.719 

0 

19 

33 

.742 

0 

20 

34 

.764 

0 

21 

35 

.787 

0 

22 

35 

.787 

180 

b) Differential Pressures 


rifice 

X( in. ) 

X/L 

^ 

8, 0-180* 

18 

.404 

0 & 

8, 90-270 

18 

.404 

90 & 


*Note: Leaky Transducer-Unusuable 
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TABLE V (cONT'd) 

FLIGHT DATA USED IN DEVELOPMENT OF LAMINAR CORRELATION 



FLIGHT DAI A USED IM DEVELOPICflT OF LAMiriAR CORRELATION 



The first step in calculatinq the boundary layers for the flight 
data is to use the measurements of static pressure along the facsimile 

cone to define theinviscid flow field. Particular data from flight tests • 

of the facsimile cone were selected to match the values of freestream Mach 
numbers listed in Table V within + 0.005. The distributions of Cp(X) were 

then used with the values of listed in Table V. to calculate static 

pressures along the surface of the BIT Cone. In the flight-test report 
by Dougherty and Fisher^^ they state that the static pressure is nearly 
constant along the length of the cone. However, the static pressure 
data we were given has significant scatter. Figures 18 and 19 show, 
respectively, a typical favorable and an adverse pressure gradient. 

In light of this apparently random scatter and the comments of Dougherty 
and Fisher, a simple straight line was fit to the data for each of the 
nine flight conditions. The associated rms errors ia Cp are listed in Table VI 
and range from 4 to 6%. Comparisons of the in-flight facsimile cone pressures 
with the theoretical distributions of Wu and Lock, Fig. 6, indicate that the 
effect of the airplane is to reduce the favorable pressure gradients 

which exist on an isolated 10-degree cone. 

The effect of orifice-induced errors was estimated using the method 
of Franklin and Wallace.^® The increase in measured static pressure, due to 
flow dipping into the 0.107 cm (0.042 in.) orifices, was estimated for two 
cases: (1) a laminar boundary layer at X = 17.8 cm and with flow conditions 
corresponding to Flight #327.0908^'^ and (2) a turbulent boundary at X = 48.3 cm 
with the same freestream conditions. In both cases, the calculated error 
in static pressure was less than 0.1%. Thus, orifice- induced errors 
appear to be negligible. 

However, the facsimile cone pressures do need to be corrected for 
differences in pitch and yaw angles. This need arises because when Mach 
numbers of flight tests of the facsimile cone and the BLT Cone are 
matched, the values of a and B do not match. A correction for these 
differences in pitch and yaw angles can be calculated via the Wu and Lock 


■^This is consistent with the Wu and Lock theory which states: Cp = 
fn(M^,, Cone angle, a, B, X/L). 

■'■^ithin the range of freestream conditions listed in Table V, this 
corresponds to the lowest Mach number for which 5 ®’'® 

measured, and thus the thinnest boundary layer ex sts for this ^ ^ 

However, in retrospect, the flow conditions for 1351 should 

have been used because a larger tw results, and orifice-induced e 

increase with 
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figure 18 
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Figure 19 


lluM)r,v. In partiailar, C^^(x) can bo calculatod nlonti tho location of 
llio I’roston-t.uho trav<’rso (il'-"0") for oach V moasurori durimi tests of tho 
laisiinilo cono and the HIT Cone. In the case of the HIT Cone, the tinie- 
averaijed value of r, diirimj the time the Prcston-tube was within tho 
laminar boundary layer, is used for tho theoretical calculation of 
(M , 1‘, X/L). The y-intercepts of the straight line curve fits of the 
facsimile cone pressures are then increased or decreased according to the 
difference between spatially-averaged (over length of laminar boundary 
layer) value of for the BLT Cone and the corresponding Cp calculated 
for the facsimile cone at the same but different r, i.e.; letting 
(Cp},:(, - mx > b, 

then (Cp)p(^ " >«x + b + Cp(V for BLT Cone) - Cp(r for FC) . (23) 

The resulting equations are tabulated in Table VI. It is relevant to 
note this procedure assumes the differences in static pressure distribution 
caused by differences in angl es-of-attack are not influenced by the F-15 
aircraft. 

The measured wall temperatures, tabulated in Table V, were obtained 

with a single flush-mounted thermocouple at x = 90.2 cm and - 180 degrees. 

Since this location is always within the turbulent boundary layer and the 

subject of this work is laminar boundary layers, it was decided to estimate 

what t hecorresponding wall temperatures should be under the various laminar 

boundary layers. This can be done by multiplying tho measured T^ by the 

ratii' of the laminar over the turbulent adiabatic wall temperature. A 

turbulenl recovery factor of r^ - 0.8825 and a laminar recovery factor of 

r 0.845 were selected because they are average values based on data from 

tests of cones, e.g. Rohsenow and Choi’*. The local edge Mach number at 

X - 80.8 cm is computed based on the corrected Cp(X) distribution ior a given 

set 01 (light data. Table VI. A laminar value for is obtained using the 

spat i<il ly-averat|etl value of Op over the length of the measured laminar Preston 

tube data. Tho measured wall temperature is then multiplied by the ratio 

(1 ) /(T ). in order to arrive at the temperatures listed in Table V. As 

^'aw'' ' awM 

may In' seen, the wall t emiiera Lures ijenerally decrease by only 2 F. 

At this point, sufficient information exists in order to proceed with 
luMiiuiary layer i om|)iita t ions for the nine fliglit conditions. This, of course 
V.as dniii' using the STAN-5 Code. The rosuUiiU) theoretical boundary layer 
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TABLE VI. Straight-Line Fits of Surface Pressure 
Distributions Measured on Facsimile 
Cone (FLT #358) 






profiles cm then he nsnrt In nonj.mctinn with I'roston-tulK prossuros to 
, I, -fine offoctivo probe heiphts. A typical trace of Prostoo-tube pressure 
is shown in fi<|. ?0, Thu scatter in the laminar portion of the liouiulary 
layer can be eliminate.1 by han.l fairinq a curve thrmuih the measured pressure. 
This has been i.1ono for all nine flight conditions. 

Using values of the faired Preston-tube pressures, at intervals of 
1.27 cm (O.G in.), a total of 87 values of were calculated. These 
are listed in Table V and are plotted in Tig. 21 as a function of 
U h/v . This figure reveals that the semi-empirical values of are 
mutually inconsistent, and values greater than two are suspected of being 

erroneous. 

However, before attempting any corrections, it is instructive to 
derive a correlation based on the unmodified data and then to calculate the 
corresponding effective unit Reynolds numbers for the 11-Ft TWT, Before 
doing this the distributions of shown in Fig. ?1, were fit with 

exponential eguations to further smooth the data and facilitate the 
corrections which will be discussed later. These equations are listed in 

Table VI 1. 

^ least-squares curve fit of a quadratic correlation to the laminar 
night data, based on the smoothed values of; (1) Cp(X), (2) 

(3) K resulted in 

Y* - 0.0915{X*)'- 0.5846 X* + 3.2259 , (24) 

for 5.8 < X* 6.4. 

A graph of this equation and the associated data are presented in Fig. 22. 
The corresponding values of are shown in Fig. 23. The rms error in 
f is 0.56:;. This reduced scatter, compared to the wind tunnel data, 
i^not only due to the use of smoothed data but also the fewer number 
oi individual data points and the smaller range of freestream conditions. 

This identifies a deficiency of the method for calculating which 
torces the measured Proston-tube pressures to be consistent with the 
theoretical solutions of the conservation equations for tho boundary layer. 
This procedure result^; in values of X* and Y* for the high values of to 
.give witi, and correlate with the rest of the data. Only the plot of 
poittts to a possible problem with some of the measured flight data. 

ihe correlation for the unmodified flight data. Kg. (24), is 
,omi>ared with the correlation for the shifted wind-tunnel data, Tq. (22), 
in 1 i.,. 24, (ho two correlations cross near X* - 5.77, which is outside 
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Figure 21 




TABLE VII 


Or<(U(WAL 

OF POOH gUAUl Y 




CO 

cn 

00 

VD 

03 

00 

LO 

CVJ 




CO 

CO 


VO 


o 


CVJ 

CVJ 

«s: 



LD 


00 

r-K 

00 

r** 

o» 

00 

p 


00 


00 

l£> 

o 

CO 

00 

00 


< 


• 

« 

* 

* 

» 

« 

* 

« 

• 

a 


r- • 



* 

CVJ 

CVJ 




H- 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

CD 



m 


1“^ 

I— 1 


1 — > 


r— 1 














c** 

CO 




cn 

VO 

<M 

00 

Uu 


CVJ 

CVJ 

o 

VO 


LT> 

CO 

o 

r— 

U_ 


cn 


CO 

CO 



00 

CO 

CO 

o oc 


o 

a\ 

CO 

lO 


CO 

o 

VO 

r— 

o 

M- 

• 

t 

« 

• 

* 

« 

« 

• 

* 

1/) u. 

0) 

00 

00 

ro 



CVJ 

o 

00 

Ov 


x: 

ro 

CO 

CO 

ro 

CO 

CO 


CO 

ro 

1—1 CO 











u. H- 

j- 

1 

1 

1 

— 1 

1 

1 

( 

1 

1 

3: 

o 










Ui CD 

M- 

H 

H 

H 

H 

H 

■ H 

‘ p 

H 

' H 

> »— * 


OC 


cc 

no 

o: 

cc 

OC 

a: 

QC 

OC LU 

c 






•x— ' 


*■>— ' 


3 ZC 

o 










o 

»r* 

OJ 

CO 



03 

CO 

LO 

00 


LU 

4J 

U) 

03 

LO 

00 

o 

VO 

r“ 



</) 00 

Ki 


CVJ 



CO 

CVJ 

CO 

o 

CTi 

UI o 

3 


r*“ 

r-* 

00 


00 

Cvl 

o 


< o- 

UJ 

o 

d 

d 

d 

d 

d 

d 

1— 

d 



1 

1 

1 

1 

1 

1 

1 

1 

1 

O' Ui 


UU 

UU 

1__l 

UU 

UU 

uu 

L^ 

uu 


c/> > 
1 1— • 


o. 

Q. 

Q- 

o. 

Q. 

a. 

CL 

CL 

Ox 

H* H- 


X 

X 

X 

X 

X 

X 

X 

X 

X 

CO o 


LU 

UJ 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

Ui U- 


o 


03 

h- 




•a- 

o 

_J u. 



00 

00 

00 

00 

00 

CO 

03 

ov 

LU 


m 

CVJ 


fO 


00 

•a- 


o 

U- 


o 

CVJ 

Lf) 

o 


o 

o 

o 

n— 

o o 


» 

• 

■ 

« 

• 

• 

• 

« 

• 

1- 


o 

o 

o 

o 

o 

o 

o 

o 

o 

CO 






















_i _j 











=3 











C/> 











UI H 











OC SZ 











UJ 











z 











o 











a. 











X 











UJ 












s 




m 

U3 

VO 

CO 

fO 




CO 

lO 

lx 

lx 

00 

00 

00 

C3V 

03 




• 

• 

4 

• 

* 

* 

• 

» 



o 

o 

o 

O 

o 

o 

o 

o 

o 



s 




00 

CVJ 

VC 

o 

lx 


fO 

o 

o 



r— 


ro 

o 

CVJ 

o 

LO 

eg 

U) 



O'. 

o 

o 


o 

03 

CO 

O 

CO 


+j 

o 


r“ 



o 

r“ 

r“ 




• 

« 

• 

• 

• 

» 

• 

■ 

* 

* 

Li. 

03 

03 

Ov 

03 

O' 


ro 

CVJ 

CO 



CVJ 

CVJ 

CVJ 

•4* 

CVJ 

CVJ 

ro 

CO 

ro 



CO 

CO 

CO 

CO 

ro 

CO 

CO 

ro 

CO 


5b 


Here = U_h/v, 

L i 



LAMINAR PRESTON-TUBE CORRELATION 
FOR UNMODIFIED FLIGHT DATA 


ORIGIi'iAL PAU;: 13 
OF POOR quality 


SCATTER OF THEORETICAL 
LAMINAR SKIN FRICTION ABOUT 
CORRELATION FOR UNMODIFIED FLIGHT DATA 


cP 

O O 


riqure 23 


6bo 







nrtorr.'j 

COMPARISON OF 

CORRELATIONS FOR SHIFTED WIND 
TUNNEL DATA AND UNMODIFIED FLIGHT DATA 


WIND TUNNEL DATA yy 

Y*=-0.0103(X^)2 +0.6653X^-0.5946 

RMS ERROR OF Of ^ =0.98% 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

// 

^ FLIGHT DATA 

j/ Y** = 0.09 15(X*)2- 0.5846 X* +3.2259 


RMS ERROR OF Cf_o=0.56% 


Range of Flight Data 

Range of W.T. Data 




the r.nc |0 of the flight data, and a maxhnon, difference ocenrs near X* = 6.06 
The wind-tunnel correlation Is above the flight data correlation over the 
range of the flight data. This Is consistent with the fact that larger 
disturbances are present within a transonic wind tunnel. The difference 
between the two correlations can be used to define an effective freestream 
unit Reynolds number. A description of this procedure Is given In the 

following section. 
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V. COMPARISON or CORIiaATIONS AND CALCULATION OF 
AN- FFFECTIVE FRFFSTRFAM UNIT REYNOLDS NUMBER 

• yojiwl if i l-U 

In analogy to the classical definition of effective freestream Reynolds 
number based on equal values of drag coefficients of a sphere, one can use 
Ei|s. {y?) and ('M) to accomplish the same objective by equating values of skin 
friction coefficient. This procedure requires the following set of conditions 
be specified at a point on the cone. 

1. Freestream conditions: M,„, q„, Re„, 

2. Inviscid pressure 

3. Preston-tube pressure 

Vf 

5. Wall temperature 

These conditions permit calculation of numbers for all the variables which 
appear in X* and Y*, except Cf. If these values are used in the correlation 
for flight data, Eq. (2A), a supposedly interference-free value of laminar 
skin friction coefficient is obtained. Now when this value of Cf and all the 
other specified variables, except Re^, are substituted into the wind-tunnel 
correlation, Eq. (22), an "effective" freestream unit Reynolds number can be 
calculated. For given freestream conditions, the effective Reynolds number 
should be independent of X. The only restriction is that the appropriate 
values of Pp, Pg, Ygff and T„ be used for a given X. However, small 
variations are found to occur with X simply due to inaccuracies in the data- 
and the correlation procedure. These small variations within a given traverse 
were eliminated by averaging over X. This has been done for all of the 19 
wind-tunnel cases, and the results are presented in Fig. 25. 

Surprisingly, the calculated effective unit Reynolds numbers are only 
about 7.5X larger than the values calculated by traditional methods. These 
I'esults are surprising because they are much smaller than those suggested by 
Reynolds numbers based on the distance to the end of transition. For example, 
data tor the 11-Ft TWT is presented in Fig. 26 which suggests the effective 
Reynolds number at M^„ 0.80 should be of the order of 35% larger than the 

normal unit Reynolds number. Also, since (Cp)rms if’ 

ll-Ft TWT, as shown in Fig. 27, and Rey is approximately proportional to 
(<’p)nns Appendix A), it appears reasonable to expect (Relief f 

peak near M = 0.70. This type of behavior was anticipated because any 
diffurences between tht? wind tunnel and flight data are assumed to be due to 
winil-tunnel-ervironment effects which were not included or modelled in the 
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data analysis, vi/., noisR and freestroain turbulence. All of this suggests 
that wind-tunnel noise alters the aver age laminar skin friction significantly 
loss than it does transition. 

However, before accepting this conclusion, we must resolve the question 
whether the semi -empirical flight values for (Fig. 21), which clearly 
point to erroneous data, have any significant effect on this conclusion. 

Thus, the question arises: Is there a rational procedure which can be 

employed to correct the Preston-tube pressures measured during the flight 
tests and thereby permit the calculation of more correct effective unit 
Reynolds numbers? The following section discusses our search for such a 
procedure. 

B. Search For A Procedure to Correct Flight Data 

If we assume Kg|r|r should be less than. two, the conclusions from Fig. 21 is 
that the data for M«. = 0.85, 0.86, 0.74 are too-h-igh. The source of these 
errors is unknown. The level of Pp for = 0.85 and 0.86 would need to be 

reduced by approximately 2.54 kPa {53 psf) in order to reduce the largest 
value of to 2. In the case of = 0.74, Pp would have to be reduced by 
2.39 kPa (50 psf) in order to match the pressure measured for the case = 
0.75, 

In addition to level of pressure, there are at least three other 
variables which must be considered. Firstly, the data is expected to follow 
the wind-tunnel pattern of decreasing Kg^-^ data with increasing M,,, , and 
secondly, orderly spacing between individual cases is expected. Thirdly, in 
any attempt to correct the Kg|:|: data, one must decide if the measured 
differences in Pp for a given trave«e are valid or not. If the errors in Pp 
are caused by 2 ero shifts in the transducer read-out equipment,, it would 
normally be correct to assume 2 ero drift was negligible during a given 
traverse w.hich required loss than a minute. However, even the differences in 
Pp for the various traverses do not agree, see Fig_2JL For example, the 
total difference in Pp between 40.6 cm X <_ 50.8 cm of Flight No. 329.1036 
(M„. * 0.74) is 0.618 kPa (12.9 psf); whereas, the corresponding difference 
for Might No. 349.1400 (M,„ - 0.75) is 0.551 kPa (11.5 psf). As shown on 
Fig. 21 and in Table V, ttiore are no significant differences in unit Reynolds 
number or freestreaiii dynamic pressure. In summary, there does not appear to 
be a concensus on either level or distribution of Preston-tube pressures! 

Ttie first attempt to correct the values of Pp, measured during flight 
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tests, is based on the idea that the flight data and the wind-tunnel data 
have a common asymptotic value of as and the Preston tube moves 

to the tip of the cone and out of the boundary layer. This requirement can be 
imposed by using the exponential equation developed from a least-squares curve 
fit of the wind-tunnel values of between = 0,6 - 0,95 (see Fig, 14) 
and the equations of Table VII, The wind-tunnel equation defines the 
asymptotic values of to be given by 0,655 (1-M^) “*^’^when U^h^y^-> <« . 
This equation is then used to define new values for the additive constants, 
which appear in Table VII, by inserting a given flight Mach number. The 
resulting modified values of Kg^f are plotted in Fig. 28, Now the Kgff are 
much lower but the order of the curves with respect to M„is still incor- 
rect, The correlation which results from use of these values of Kgff is 

Y* = - 0,00016 (X*)2 + 0.5054X* - 0.0194 . (25) 

This equation is plotted in Fig. 29 and is compared with the shifted 
wind-tunnel correlation, Eq. (22). The associated rms error in is now 
0.15%. This extremely low error in Cf is caused by the close proximity to the 
wall. In this region, total pressure within the laminar boundary layers is 
relatively insensitive to changes in Kgff, see Fig. 13, and the corrected 
values of Pp are approximately equal to P„. This results in nearly constant 
values of Pp and is considered to be invalid. However, the corresponding 
effective Reynolds numbers vary nearly monotonically from a maximum of 18% 
above the normal Reynolds number for =0.30, Re^ = 13.1 x 10'’ down to a 
minimum of 4% at M^^, = 0.90, Re^ = 9.8 x 10^. A graph of these Reynolds 
numbers is presented in Fig. 30. 

In light of the above inconclusive results, a second correction pro- 
cedure was investigated. This second procedure is based on the assumption 
that the effects of freestream turbulence and noise on the measured laminar 
Preston-tube pressures are a minimum at the beginning of boundary layer 
transition. This, in effect, assumes the flight and wind tunnel values of 
Kgff are equal at X^- for a given R^ and Mo,. This then leads to use of the 
asymptotic curve fit of the shifted wind-tunnel data. Fig. 14, to calculate 
Keff(Xt), The theoretical values of R^ - at X^ are substituted into 

Kprr = 2.865e‘®-®2” R + 0.655 (l-M^) i (26) 
for the nine different flights. The differences between these new values of 
Kgff (Xf) and the original flight values defines an increment which was added 
(or subtracted) from the complete set of Kgff's for a given flight. The 
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resulting correlation is 

Y* = 0.1247(X*)^ - 0.9210 X* + 4.0548, 
for 5.6 < X* < 6.0. 

This equation is compared in Fig. 31 with the preferred wind-tunnel 
correlation, Eq, (22). Unfortunately, the two correlations cross. 

This results in effective Reynolds numbers varying between 65^ above 
and 6% below the normal Reynolds number. Thus, this procedure is 
not helpful . 

A third correction procedure focuses on using the unmodified 

values of for Flight No. 349.1400 = 0.75) as a reference for 

the rest of the flight data. This case is attractive because it appears 

to have the most realistic values of i.e., 1.7 < <1.8. In 

order to apply this assumption, it is necessary to determine-bow the 

rest of the data should be distributed about this reference 

case. This can also be done using the curve fit of the shifted wind- 

tunnel data, Eq. (26). Corrections for the three cases in the upper 

left of Fig. 21 were calculated by using the exponential curve fits 

for each case (Table VII) to define a corrected at a value of 

corresponding to the midpoint of the reference case, = 0.75. 

The corrections for the five cases on the right of Fig. 21 were 

calculated by using the exponential curve fit of the reference case 

to extrapolate to an R^ defined by the midpoint of each of the five 

cases. The reference Mach number, 0.75, and the five different values 

of R^ were used-Jn Eq. (26) to calculate reference values of The 

five different values of R and M were then used in Eq. (26) to define 

00 ' ' * 

"correct" differences between these data and the reference case. Each 
of the calculated differences in were then used to shift the center 
of the data for the corresponding eight flights. This procedure 
provides a more realistic spacing of the eight flights about the 
selected reference case. In addition, this time it was assumed the 
measured differences in Pp for a given traverse are valid and should be 
maintained. This requirement leads to a particular distribution of 
about the shifted midpoints. The results are illustrated in Fig. 32. 

The rather weird distributions of are caused by the fact that is a 
function of both x and y and the derivative dP^/dx changes with shifts in 
distance from the wall, y. The associated correlation is 

Y* = 0.0598(X*)' - 0.1777 X* + 1.^28, (28) 

for 5.64 • X* < 6.09 . 
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Thif. equation Is shown in Fig. 33 along with tho individual data Points. 

A c:nmpar ison botwoon Fq. and tho preferred wind-tunnel correlation is 
Micido in Mg. 34. 

The small scatter of theoretical skin friction coefficient about values 
computed from Eq. (20) is presented in Fig. 35. The rms error in skin 
friction coefficient is 0.37% which is less than the 0.97% scatter found for 
the shifted wind-tunnel data, Fig. 16, It is revelant to here note that the 
scatter in generally decreases as Y^ff decreases. This is illustrated by 
comparing the values of K^ff shown in Fig. 28 with the higher values of Fig. 

32 and noting the lower rms error in of 0.15% as indicated_on Fig. 29. 
Thus, one can legitimately question whether the flight data have been shifted 
down too much toward the wall. 

Before addressing this question, we need to first present the values of 
effective unit Reynolds numbers which result from using the two correlations 
shown in Fig. 34. Values for the nineteen wind-tunnel conditions are shown in 
Fig. 36. The maximum effective Reynolds numbers are approximately 6.5% larger 
than the corresponding standard wind-tunnel values. These are about 1%-less 
than those obtained for the unmodified flight data, Fig. 25. The distribution 
lias only a vague resemblance to the measured noise curve shown in Fig. 27. 

The two lowest values of effective Reynolds number (~1.04 Re^p) occur for 

- 0.3, Re„^ = 13.1 X lO" and = 0.9, Re^ = 9.8 X lOS Data. for both of 
these two cases depart from the majority of the data shown in Fig. 14. The 
low Mach number case appears to be too high and has values of shifted 
toward the flight data. The corresponding values of X* lie between 6.26 - 
6.28. As indicated in Fig. 34, the difference between the two correlations is 
less in this range of X* and this results in a smaller value of effective 
Reynolds number. The high Mach number case appears to have a slope dKg^r^/dR 
that is conspicuously large. However, the primary reason Re„,^gff is small for 
this case is that X* lies in the range 5.73 - 5.80, This corresponds to the 
far loft: of Fig. 34 where a smaller difference between the two correlations 
also exists. Finally, we face the questions what happens to the effective 
Reynolds numbers if the flight values of Kgff are shifted either up or down. 
This was investigated by shifting the data of Fig. 32 according to the 
fnl1()winq procedure. The average R .j. and the corresponding was calculated 
for the I'eference case (M_„= 0.75) and likewise for the other flight 
conditions. The data were then shifted by arbitrarily adding a constant 
increnuMit to each of the average Kg^r^'s. These new distances from the wall 
were then used tolocate a new value of Pp within the corresponding STAN-5 
Ijoundary layer at that streamwise station. The increment in Preston- 
;■> 
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tube pressure, for a given flight condition, was then added or subtracted from 
all of the Pp for. that traverse. These new values of Pp were then used in 
conjunction with the theoretical profiles to define new distributions of Kgff 
and the associated values of X* and Y*. A new least-squares curve fit of 
these data provides a 

new correlation equation which can then be used with the wind-tunnel 
correlation to define a new set of effective unit Reynolds numbers. 

The reference value of =1.72 was increased up to 2.12 and was 
decreased down to 1.42. The resulting distribution of Re^^gff is presented in 
Fig. 37 for = 0.80. This clearly shows that the calculated effective 
Reynolds number have a maximum of approximately 1.065 Re,„. Thus, we may 
conclude that the procedure selected for correcting the flight values of 
reduces the maximum values of Re,„^gff by about 1% of Re„,. ^ 
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CONCLUSIONS AND RLCOMMCNDATIONS 


Analyses of Preston-tube data, obtained with the AEDC Transition Cone 
in both the 11 -ft TWT and flight tests, have revealed that these data have 
errors. The errors in the 11-ft TWT data are relatively minor and are 
comparatively easy to correct. A satisfactory correlation between the 
corrected wind-tunnel Preston-tube pressures and theoretical laminar 
skin friction has been developed. The standard deviation between the cor- 
related skin friction coefficients and the corresponding theoretical values 

is 0.98%. 

The errors in the flight data were found to be more severe and a 
great deal of effort was expended in the search for a rational correction 
procedure. Three distinct correction procedures were investigated, and 
correlations have been developed based on each of these modified sets of 
flight data. The preferred correction procedure forces the flight data 
to exhibit some of the orderly characteristics of the wind-tunnel data. 

The corresponding preferred Preston-tube/skin friction correlation 
exhibits an rms error in skin friction of only 0.37%. 

The wind-tunnel and flight correlations have been successfully used 
to define "effective" freestream unit Reynolds numbers for the 11-ft TWT. 

Based on the preferred rearrangement of the flight data, the maximum 
effective values of Re^^ are approximately 6.5% higher than the normal 
tunnel values. This compares with 7.5% which results when using the 
unmodified flight data. The maximmii effective Reynolds numbers occur 
for Mach numbers between 0.60 and 0.80 and for nomal tunnel unit Reynolds 
numbers of 9.8, 13.1, and 16.4 million per meter. 

The distribution of (Re,„)^^^ with Mach number only vaguely resembles 
the distribution of total noise intensity as measured in the 11-ft TWT 
with microphones on the AEDC Cone. Also, the values of {Re,^)eff 
for a given noise intensity when compared to the noise/boundary-layer- 
transition correlation of Doughterty and Fisher . Thus, we are led to 
conclude that wind-tunnel noise effects the aye^ laminar skin friction 
much less than it effects boundary layer transition. 

It is reconinended that the analyses used herein to develop Preston-tube/ 
skin-friction correlations be applied next to the analysis of turbulent 
boundary layers and finally to the transition region which is known from 
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Mk' work of Dougherty, et to be sensitive to tunnel noise. These 

new correlations can be used to define additional values of (Re^^leff which 
can be compared to values reported here, 

Boundary Layer Transition Data 

Date for the onset, extent, and end of boundary layer transition are 

summarized and presented in a systematic form in Appendix A. Boundary-layer- 

transition data for the 11-ft TVIT appear to be correlated by the product 

Re (C In the range 0.60 < M < 0.95, this product is near 4 and 

T p niis 

decreases rapidly as Mach number decreases below 0.60. 

Effects of Nose Radius 

The effects of variations in nose radius on the transonic flow about 
a cone are discussed in Appendix B. The radius of a hemispherical nose has 
been found to have an imporlaat and dominant effect on boundary-layer 
transition and flow separation at supercritical speeds. Increasing nose 
radius allows the growth of a thicker boundary layer prior to encountering 
the transonic shock near the juncture of the nose with the cone. -This 
promotes both earlier transition and laminar flow separation. An analytical 
procedure has been developed to model this phenomena, and a comparative 
study of the effects of nose radius, Mach number, Reynolds number, and heat 
transfer can now be perfomed. 
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Al'PENDIX A 

SUMMARY OF SUBSONIC BOUNDARY LAYER 
TRANSITION DATA 


W ind Tunnel Data 

The surface distances from the nose of the cone to the onset (X^) 
and end-of-transtion (Xj) are tabulated in Table A-1 for twenty one 
different subsonic, wind-tunnel conditions. The test conditions include 
Mach numbers from 0.3 to 0.95 and free-stream unit Reynolds numbers of 
3, 4 and 5 million per foot. The extent of the transition zone (AX), the 
related ratio Xj/X^, and the root-mean-square of the fluctuating pressure 
coefficient [(C )y.n, 5 ] are also included in Table A-1. (Cp )^,„5 is based on 

microphone data^taken with a 1/4 in. diameter microphone mounted flush with 
the cone surface at.a. distance of 18 in. aft of the nose. Data over the 
frequency spectrum between 0.2 and 30 kHz were used to define • 

Transition Reynolds numbers, based on the product of free-stream unit 
Reynolds number and the distance to the end-of-transition, are included be- 
cause this parameter was initially used by Dougherty’ to correlate with 

(C 

^ ’TH^'tfleld and Dougherty^’ have stated that a wide variety of previous 
test data indicate the ratio of X.j./X^ is approximately 2 and is nearly 
independent of flow conditions. However, the values of this parameter in 
Table A-1 range from 1.10 to 1.23. There appears to be a slight tendency 
for Xy/X^ to increase with Mach number. As is shown later in Table A-2, 
the flight values of X.^/X^ also fall within this same range. Thus, the 
value of two is not characteristic of these tests. 

The effects of M^ and free-stream unit Reynolds number on X^ and X^ 
are shown, respectively, in Figs. A-1 and A-2. The primary results to note 
are that the distances to onset and end-of-transition (1) increase with 
decreasing Re^^, but vary only slightly with M^ and (2) are smaller, i.e., 
occur closer to the nose, when M^^ = 0.70. This corresponds to the Mach nuin 
ber at which (C reaches a maximum in the Ames 11 -ft TWT, e.g., see 
Dougherty and Fisher.” In Ref. 29 , Dougherty presented a correlation of 
wind-tunnel data which indicated Re.^. = 3.7 x 10® (Cp)“^ 5 » within + 20%. 
Thus, Rej(Cp);^s plotted in Fig. A-3 as a function of M^^ in order to 
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PRESSURE COEFFICIENT: 11-FT TWT DATA 








ascertain if this parainotor Is Indood Independent of Mach number, for 
between 0.60 and 0.95, an average value of approximately 4 x 10^’ appears to 
be a good correlation. In this Mach number range tunnel -generated noise 
reaches the highest Intensity levels. As M^, decreases below 0.6, the noise 
level also decreases rapidly. This rapid decay In noise Is reflected in 
the decreasing values of Rej(Cp)^^g as Mach number decreases. The data 
suggests that this parameter is a unique function-of M^ for the AEDC Tran- 
sition Cone In the 11 -FT TWT. However, there Is Insufficient data to support 
a definite conclusion. 

The effects of M^ and Re^^ on the length of the transition zone are pre- 
sented in Fig. A-4. The exteni.-of-transit1on appears to have a local mini- 
mum near M^ = 0.70 only for the case Re^^ = 3 x 10®. The significance of 
this observation in unknown at-±his time. 

Before proceeding to discuss the subsonic flight data, it is relevant 
to note that wind-tunnel data were used to estimate corrections of and X-j. 
for flight cases which had nonzero pitch and/or yaw angles. Figures A-5 
and A-6 present the data for.. AX as a function of r, the total angle-of-attack; 
the corresponding vaues of X^ and X^ were used to correct the subsonic flight 
data. Figure A-5 expresses the data as a nondimensional ratio AXj./AXj. ^ q vs. 
r, and Figure A-6 recast the data in terms of (AXp - AXp _ wind- 

tunnel data were obtained from tests of the 10-deg cone in four different 
tunnels. In all cases, the pitch angle (a) was varied with yaw angle (p) 
equal to zero, and then p was varied with a = 0. 

The flight values of a and 3 are known from measurements with the cali- 
brated flow angularity probe (yawmeter). The following equations were used 
to calculate values of r and i}) (azimuthal angle of the Preston tube, measured 
positively in the clockwise direction from the windward element of the cone, 
looking forward along the axis). 

r = («^ + 3^)^ 

= tan-‘ - 11 

These values were then used to interpolate wind-tunnel values of X^/(X^)j. , 
and X|/(Xy)j, _ ^ as a function of r and 4>. The measured flight values of 
X^ and X-|. were corrected for nonzero r by dividing by the ratios which 
corresponded to an M_, nearest to the flight Mach number. 


A-6 


ORICJt'JAl. PAGE IS 
OF POOR QUALITY 

EFFECTS OF Mod AND Reft ON EXTENT OF 
TRANSITION: 11 -FT TWT DATA 





Qt- VO'jU quality 


EFFECTS OF YAW AND PITCH ANGLES 
ON EXTENT OF TRANSITION: RATIO OF AXpTO AXp,© 


/ ° 


f 

4 





^ WINDWARD^ 




SOLID SYMBOLS DENOTE DATA 
AT 90** FROM WINDWARD 



Mco 

Reft 

N.L. 

N.A. 

AEDC 

A 

04 

3.0M 

16TDT - - 



□ 

0.6 

3.0M 

— 

14TWT 

— 

V 

0.8 

3.0M 

— 

— 

16T 

O 

0.9 

3.0M 

— 

11TWT 

— 

O: 

0.9 

3.8M 

— 

14TWT 

— 

o 

0.9S 3.0M 

— 

11TWT 

— 


Figure A-5 
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Figure A-6 






The inconsistency of the data shown in Figs. A-5 and A-6 suggests this 
procedure becomes progressively more inaccurate as r increases and as the 
difference between flight and wind-tunnel Mach number increases. Fortunately, 
the subsonic values of r v/ere of the order 1/2 deg or less, see Table A-2. 

Only the data from the Ames 11-Ft and 14-Ft tunnels at = 0.90 are comparable. 
When r = 0.5 deg, Fig. A-5 shows a difference in. AX on the leeward side of 
0.28 (AX)p _ Q or a 28% discrepancy; whereas on the windward side, tiie dif- 
ference isl4%. An estimate of the errors for other Mach numbers is not 
possible. Because of this uncertainty in the procedure for correcting flight 
transition data for nonzero Land the inability of the SIAtl5 program to model 
asyiranetric boundary layers, only flight cases for which [rj < 0.11 deg were 
used in the-development of the Preston-tube/skin-friction correlation. 


Flight Data 

The subsonic flight data are summarized in TaJale A-2. As discussed pre- 
viously, the last four digits in the flight number designate time of day dur- 
ing which the data were taken. For example, 349. .1347 denotes flight number 
349, and the data were obtained at 13:47 hours. In addition to nonzero values 
of r, the "corrected" data for X^, X.j. and AX have also been modified to 
account for non-adiabatic wall temperatures via the following equation. 

^’^t^Tj^T _ ^^T^Tj^T 

^ aw a 

This is an empirical equation which is based on a curve fit of flight data, 
see Ref. 8. 

Figures A-7 and A-8 present the distance to transition onset for the 
uncorrected and corrected data as a function of M^ and Re^^. Figures A-9 
and A-10 present similar plots of distance to end-of-transition for the un- 
corrected and corrected data. In both cases, the corrected values of X^ and 
Xj show a considerable improvement in defining a discernible dependence on 
M^ and Re^^. The general trend for both X^ and X^ is that transition occurs 
earlier {i.e., X's decrease) with increasing Re^^ and decreasing M^. 

In Fig. A-11 transition Reynolds numbers Re^, based on the product of 
Re ft and Xj, are shov/n as a function of M^. It is relevant to here note 
that the difference between and for the subsonic flight conditions 

is less than 1%. Since the data is not collasped into a single curve, it 
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Figure 




appears that other variables, such as are involved. Unfortunately, 

(Cp)^ii^g was not recorded during the majority of the subsonic flights. 

The extent of the transition zone is plotted in i'igs. A-12 and A-ld. 

In this case, the corrected values in Fig. A-13 show no discernible reduc- 
tion in scatter, and the correlation remains unsatisfactory. 
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AI’I'KNDIX B 


lIl-fECTS OF SPHERICAL NOSE BLDNTNESS ON 
BOUNDARY- LAYER TRANSITION AT SUPERCRITICAL SPEEDS 

In Ref. 8, Dougherty and Fisher declare that the AEDC-BLT Cone has an 
apex with an "equivalent" diameter of 0. 004-in. (the word equivalent is 
apparently used to denote the fact that the nose-tip is not exactly spheri- 
cal.) Since it is known that hemisphere-cylinders generate rather strong 
shock waves at high subsonic Mach numbers, e.g., Hsieh®°, a natural 
question is; how does nose bluntness affect boundary-layer transition with 
in the transonic regime? 

In order to answer this question via analysis, it is necessary to have 
at least two computational tools. One is a computer code which will calcu- 
late both the inviscid flow with shock waves and the associated profiles of 
the boundary layer. The second required tool is a procedure to estimate 
when and where boundary layer transition occurs. The computer code developed 
by Nietubicz, et al.'^ was identified as a state-of-the-art method for solv- 
ing the Navier-Stokes Equations about axisymmetric bodies and at transonic 
speeds. The second needed tool was located by contacting Paul Granville at 
NSRDC. He was contacted because he had previously published papers dealing 
with boundary-layer transition on axisymmetric bodies. Although none of his 
work was directly applicable to the transonic cone question, he did supply 
us with a copy of a very recent paper by Wazzan, et al.^^ 

The important results of this paper is the following equation for esti- 
mating the onset of boundary-layer transition. 

logio [Rg(e'’)] = -40.4557 + 64.8066H - 26.7538H^ + 3.3819H^ . (B-1) 

for 2.1 < H < 2.8. The authors claim that this equation correlates the well- 
known "e’ method", see Ref. 33. Thus, the use of Eq. (B-1) avoids the need 
for lengthy stability calculations. The authors state that the method is ap- 
plicable to incompressible flows which: 

1. do not vary too much from local similarity, 

2. have small surface roughness and/or vibration, 

3. low freostream turbulence, and 

4. small heating rates, i.e., T^ - T , '' 23'’C. 
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For the purpose of conducting a comparative study of the effect of noso- 
bluntness on boundary-layer transition on the conn, it was decided to apply 
Eq. (B-1) to the compressible boundary layer by evaluating the kinematic 
viscosity that appears in at the reference temperature of Sommer and Short, 
Eq. (11). 

Next, Dr. Nietubicz was approached about using his program for this type 
of analysis. He agreed-that this was a problem that no one has studied and 
thought his code should Jie able to do the job. Unfortunately, his Navier- 
Stokes solver did not print. out the conventional, boundary- layer parameters. 

He stated that he didn't have the time to add a subroutine to accomplish 
this and suggested the first author visit BRL to perform this task. A three- 
day trip was made and with some additional debugging via the mail, this sub- 
routine was successfully added to the BRL program. 

The original agreement with Nietubicz was that one or two cases would 
be run in order to check for satisfactory execution, and at that point. Dr. 

Pulliam at NASA Ames wouIcUcun the code for a set of freestream conditions 
and different nose radii. An initial case was run during the check-out pro- 
cess for a spherical nose radius of 0.05 in. and with M = 0.95, Re,^ = 3 x 10®, 

“ft 

= 450“R, and = 530^R. The results are summarized in Figures B-1 thru 
B-3. Figure B-1 shows a portion of the 40 x 80 set of grid points that were 
used in the finite difference solution of the equations. Figures B-2 shows 
the distribution of pressure coefficient about the nose. Here, the axial 
distance X is normalized with respect to a diameter of 1.4 in., and in order 
to focus on the nose region, only the first 8 in. of the cone was modelled. 
Figure B-3 presents the corresponding Mach number contours. The central in- 
terior Mach contour encloses the region with local Mach numbers ^ 1.10, and 
Mach number decreases 0.01 with each larger contour out to a Mach number of 
0.95. The pressure distribution and the Mach contours indicate a rather 
smooth recompression as the flow passes from the nose cap onto the cone, and 
furthermore, the flow remains attached. As freestream Mach number is in- 
creased, it is anticipated that a shock v/ave would form and the boundary layer 
would thicken and possibly separate as M^ approaches one. However, the above 
calculation assumes a laminar boundary layer. If transition to a turbulent 
boundary layer occurs, separation would be either delayed or, possibly, would 
not occur at all. The boundary-layer subroutine was not ready at the time 
this case was computed; thus, no information is presently available concerning 
transition. 
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A second case was computed after the boundary-layer subroutine became 
operational. This case consists of a nose radius of 0.1 in., = 0.95, 

Re = 3 X lOS = 490“R, and = 530"R. The pressure coefficients for 
this case are shown in Fig. B-4. The recompression followed by a second 
expansion is characteristic of a separated-flow region. The associated Mach 
contours are shown in Fig. B-5, and the coresponding streamlines are pre- 
sented in Fig. B-6. This figure vividly illustrates the pocket of recircula- 
ting flow. A typical velocity profile within this separated region is pre- 
sented in Fig. B-7, and a velocity profile slightly downstream of the separated 
zone is shown in Fig. B-8. The boundary-layer thickness in a separated region 
is defined to be measured from the point at which the velocity becomes posi- 
tive. 

The boundary-layer properties upto the beginning of separation are 
summarized for this case in Table B-1 . The left and right-hand sides of 
Eq. (B-1) are also include in the last two columns of this table. There are 
a number of observations that should be noted from Table B-2. Firstly, the 
boundary-layer thickness oscillates but the displacement and momentum thick- 
ness (integral properties) are monotonic. The oscillations in 6 could pro- 
bably be eliminated by placing more grid points within this region of large 
and rapidly varying pressure gradients. Finally, equality between the left 
and right-hand sides of Eq. (B-1) first occurs between X/D values of 0.790 
and 0.794. Unfortunately, the corresponding values of H are greater^than 
2.8 and thus violate the restriction on H imposed by Wazzan, et al. The 
possible error this causes in the predicted location of transition-onset is 
presently unknown. However, the pressure gradient is small in this region, 
and the assumption of local similarity is probably valid. 

Since transition onset occurs near X/D = 0.794, the flow field needs to 
be recalculated using this information. It is expected that the extent of 
the separated region will diminish considerably if not completely disappear. 
Unfortunately, the subject code does not presently have the capability to al- 
low boundary- layer transition, i.e., the flow is either entirely laminar or 
entirely turbulent. This capability is needed in order to arrive at a satis- 
factory model of the complete flow. 

In summary, the nose radius determines the wetted length and thereby the 
boundary thickness at the shock. A larger nose radius not only pronotes 
boundary-layer transition at transonic speeds but can also lead to separation 
of a laminar boundary layer. As long as one is not too concerned about precise 
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(ihsniij1;o vnliin'i, the iiiodc'l d i<:cijr>r,r!rl herein apiir-nrs to he oppropr into lot’ a 
c ompa rative study of the effects of nose-bluntnoss at transonic speeds, II 
and when the parametric study Is completed, the results will be reported 
elsewhere. 
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